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STRAIN FRAME AND THE VELOCITIES OF 

POINTS 
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Abstract 

We obtain an exact strain consistency equation for full, elastic and 
plastic strain characteristics that have a clear physical meaning and are 
naturally related to stresses. The dynamic equations are represented in a 
form that does not use the objective stress rate. This obviates a number 
of essential difficulties in the theory of finite elastoplastic strains. 



1 Introduction 

The present paper gives a detailed justification of the results obtained in [1]. To 
characterize them, let us very roughly outline the main difhculties encountered 
in mathematical elastoplasticity when passing from infinitesimal to finite strains 
and the methods used so far to overcome these difficulties. A detailed analysis 
can be found in the critical survey [2] of the state of the art in the theory by 
1990 and in the monograph [3], which summarizes what has been done in the 
last fifty years. 

Since the above-mentioned difficulties already arise when dealing with ho- 
mogeneous isotropic elastic-perfectly plastic materials, we momentarily restrict 
ourselves to such materials. 

In the Eulcrian variables, elastoplastic deformation processes are described 
in terms of mass point velocities, material density, selected measures (charac- 
teristics) of full, elastic, and plastic strains, and stresses. With the advent in 
1930 of the Prandtl-Reufi equation (see [4]), the closed system of elastoplastic- 
ity equations for infinitesimal strains was obtained. We single out the following 
three of these equations (the remaining three equations remain unchanged under 
the transition to finite strains): 



A. The strain consistency equation (a tensor equation) relates full, elastic, 
and plastic strain characteristics. 

B. The stress response function (a tensor equation) expresses the stress via 
an elastic strain characteristic. 

C. The flow law (a tensor equation) relates the plastic strain to the stresses. 
Under the transition to finite strains, the strain consistency equation (A) is 

not true in general ([3], 7.2) and the system is no longer closed. Moreover, not 
only measures but also the very notions of elastic and plastic strain have no 
precise physical meaning, and there is controversy (as yet unsettled) as to how 
exactly they should be introduced ([2], 4A). 

In this connection, there has been quite a few general theories (some of which 
are presented in [5-11]), a fact that alone shows how unsatisfactory the situation 
is. All these theories can be divided into four groups. 

1. Elastic strain is not considered at all. The plastic strain is introduced by 
a constitutive equation of type (C) via a measure that does not have a precise 
physical meaning. Two variables (plastic and full strain) are introduced in the 
constitutive equation (B) instead of elastic strain, and the system again becomes 
closed (e.g., see [5]). 

D. The fact that the strain measure is artificial casts doubt on the existence 
(and all the more, the simplicity) of the constitutive equations (B) and/or (C). 

2. The plastic strain is not considered, and so equation (C) disappears. 
Instead of (A), using various analogies, one introduces a constitutive differential 
equation for the elastic strain, which closes the system. It is this new equation 
that is doubtful (e.g., see [6, 7]). 

3. Both elastic and plastic strains are introduced. One of them is determined 
via the other (using the full strain) but has no precise physical meaning. 

In this case, the very definition already gives the strain consistency equation 
instead of (A) . The system remains closed but has the drawback described above 
in (D) (e.g., see [8, 9]). 

4- For the strain consistency equation (A), one uses the multiplicative de- 
composition of the full strain gradient in terms of the elastic and plastic strain 
gradients or the related additive decomposition of the full strain rate in the 
elastic and plastic strain rates [10]. 

Serious objections to the first relation can be foimd in [2], 4A. The second 
relation is not true in general and involves the objective stress rate, resulting in 
heavy inconsistencies ([2], 4F; [3], 7.2, 7.3). 

In the present paper, elastic strain is characterized by the elastic strain frame 
a, which is more informative than the Cauchy-Lagrange tensor. In the deformed 
medium, the elastic strain frame field generates the Riemannian metric of the 
"natural, imstressed state." Under sufficiently wide assumptions, this metric 
permits one to judge about the natural measure of deformed material objects 
(i.e., the measure that they would have after unloading), without resorting to 
actual unloading. 

We introduce the plastic strain tensor, which has the same physical meaning 
in terms of the natural measure as does the elastic strain tensor for elastic media 
in terms of the geometric measure. 
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We obtain the strain consistency equation in frame form, providing an exact 
relationship between the full, elastic, and plastic strain characteristics. 

We introduce the relative stress tensor (with respect to the Riemannian 
metric generated by an arbitrary nondegenerate frame field), whose restrictions 
are the natural stress tensor and the conventional Cauchy stress tensor and 
Piola-Kirchhoff second stress tensor. 

For media homogeneous and isotropic in the natural state, we suggest a 
plastic constitutive law specified by a scalar function p of the principal stresses 
and possibly by parameters determined by the strain history. The presence or 
absence of yield surfaces is determined by the properties of p. 

Under the assumptions accepted in the paper, we obtain a system of quasi- 
linear first-order partial differential equations in Eulerian variables for v and a, 
describing the elastoplastic deformation process in the entire space-time domain 
occupied by the medium. All equations are given in the form solved for the total 
time derivative, which makes the system suitable for numerical solution^ 

2 Some notation and conventions 

We assume that the real space in which the continuous medium moves is the 
three-dimensional Euclidean affine space. We denote it by E and the associated 
vector space by E. Let (•, •) and | • | be the scalar product and the corresponding 
norm in E. 

A frame consisting of vectors G E, i = 1,2,3, will be denoted by d 
and written as a column matrix, so that d := (01,02,03)"^, where T stands 
for transposition. (For convenience, the entries of a row matrix are separated 
by commas.) In the presence of an ambiguity, we always take a nondegenerate 
right frame. The symbols := and =: stand for " is by definition equal to", where 
the object to be defined is written next to the colon. 

Let us define the norm of a frame in E by the formula ||d|| := (laiP + |£i2p + 
losl^)"'^^^ and the difference of frames by the formula d — b := (ai — bi,a2 — 
b2,a,3 — bs)'^ . The coefficients (contravariant coordinates) in the expansion of 
a vector x G E in the vectors of a frame a will be denoted by a;°, or, more 
briefly, by xf; x°- := := (a;°, Xj, 0:3)'^. Thus x = xfai. Here and in what 
follows, summation from 1 to 3 is assumed over a repeated nonunderlined index 
in a monomial. 

The scalar product defined on E can be represented in coordinate form in 

^Note that the elastic strain frame was introduced for elastic media in I. A. Solomeshch 
and V. Sh. Khalilov, Nonlinear Elasticity Equations in Strain Frame Components [in Russian], 
No. 5089-B90, VINITI, Moscow, 1990 and for elastoplastic media in I. A. Solomeshch and 
M. A. Solomeshch, Dynamic Elastoplasticity Equations for the Elastic Strain Frame and the 
Velocities of Points [in Russian], No. 214-B95, VINITI, Moscow, 1995. The latter paper 
indicates the possible connection of the elastic strain frame in polycrystal materials with the 
strain of the crystal lattice. In the earlier-mentioned papers [6, 7, 9], frames actually coinciding 
with the elastic strain frame are introduced under various names (on the average, the trivector 
of microstructure variables) to characterize the elastic strain in such materials. 
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any frame e orthonormal with respect to this product: 

{oc,0)=atf3f Va,f3eE. (2.1) 

Conversely, for an arbitrary nondegenerate frame a one can construct a scalar 
product (generally different from (•,■)) in E by setting 

(a,/3)^ :=</3f; (2.2) 

the frame d is orthonormal with respect to this scalar product. The Euclidean 
affine space obtained by the replacement of the scalar product (•, •) by (•, •)a will 
be denoted by Obviously, Eg = i? if the frame e is orthonormal in E. 
The modulus of a vector a lE E in Ea is equal to 

\a\a := {cx,a)y^. (2.3) 

The volume of the parallelepiped spanned by a frame b in the space E will be 
denoted by and the volume of the same parallelepiped in Eg, will be denoted 

by \b\a. 

Nondegenerate frames d and d' generating the same scalar product, i.e., such 
that {a,(3)a = {cx,/3)a' for any cx, /3 G E, will be called equivalent. 

For frames d and b, by we denote the matrix whose jth column is a^. 
Matrices are sometimes indicated by the symbol " . For column matrices c := 
(ci, C2, cs)^ and d := (di, ^2, ds)-^ with numerical entries, we write 

{c,d):=Cidi, |c|:=(c,c)V2 (2.4) 

Let a continuous medium occupying an open spatial domain V c E a.t the 

initial time to fill a domain Vt at the final time t. The spatial position of 
a particle M at time to will be denoted by x (the Lagrangian variable), and 
the position of the same particle at time t will be denoted by x (the Eulerian 
variable) . 

The function x = x{x) specifying the correspondence between the initial and 
final positions of a mass point, as well as the laws of motion considered in what 
follows, is assumed to be one-to-one and have an invertible derivative x'^. 

In general, we adhere to the following principle whenever possible: objects 
characterizing the medium at the terminal time are denoted by italic Latin 
letters, and the corresponding objects at the initial time are denoted by the 
same letters typeset in Roman; any function of a particle expressed via a; or x is 
denoted by same letter with the argument indicated where necessary; notation 
introduced for some object is later used for other objects of the same nature. 
Unless explicitly specified otherwise, we assume that all functions occurring in 
the text (including the derivatives) are continuous, all curves and surfaces are 
smooth, and all sets over which integration is carried out are measurable and 
closed. Whenever the metric is not indicated explicitly, the metric of the main 
space E is meant. 

In the following, to reveal elastic and plastic strains of the medium, we 
need to compare the measure of material objects in Riemannian spaces that are 
different in general. The corresponding mathematical technique is given in the 
following section. 
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3 The change in the measure of geometric ob- 
jects under a mapping of Riemannian afRne 
spaces 

3.1. The Riemannian metric generated by a frame field. Let V be an open 
domain in an affinc spacc^ E. (Momentarily, it is not important whether E is 
Euchdean.) A Riemannian metric on V is usuaUy defined with the help of a 
metric form 

9ij{x{0)d^id^i (3.1) 

in some curvilinear coordinate system x = .t(^). Essentially, formula (3.1) 
defines a scalar product at any point x V (more precisely, on the tangent 
space at x, i.e., in our case, on E). Therefore, one can define a Riemannian 
metric by directly specifying a scalar product at the points of V. 

For our purposes, it is convenient to specify each of these products by spec- 
ifying a frame orthonormal with respect to this product. 

Thus we adopt the following scheme. Let a nondegenerate frame field a{x) 
be given on V. For the scalar product at a point x we take the product (•, •)a(x)j 
with respect to which a{x) is orthonormal (see (2.2)). The Riemannian space 
thus defined is said to be generated by the frame field a{x) and will be denoted 
by Va- The metric in it will be called the a-metric. 

The relationship between the coeflicients of the metric form and a{x) is as 
follows: 

where 

dx ( dx dx dx ^ 



Note that k{^) is the local frame of the curvilinear coordinate system at the 
point a;(^). 

The conventional formulas of Riemannian geometry in terms of the frame 
Q.{x) acquire the following form. The cosine of the angle between vectors x and 
y at a point x is 

/' — -^ {•^Ty)&{x) /o o\ 

cos{x,y)^^^) := — — ; (3.3) 

\x\&(x) ■ \y\&(x) 

The length of a curve I: x = x{t), t e [tq, n], is 

\l\&-= [ \x' {T)\a(xir)) dr; (3.4) 

J To 

The area of a surface s: x = x{u), u = {ui,U2)^ G a C R^, is 



S\& - \/\^'ui\l(x) ' \^'u2\&(x) ~ (^'ui'K2)&(x)^'^' (3-5) 
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The volume of a domain v. x = x{^), ^ = (^i, ^2, ^3)^ e w C il^, is 



/ \m\ai.md^- (3-6) 



If 6i{x) is a constant do, then V^o is the Euchdean space with scalar product 

3.2. The dilatation coefficients under a change of the Riemannian metric. 
Let a{x) and 0{x) be nondegcnerate frame fields defined on V, and let 
and be the Riemannian spaces generated by these frame fields. Let us find 
the dilatation coefficients for the measures of one- to three-dimensional objects 
under the transition from their measurement in the d-metric to the measurement 
in the /3-metric. 

The length dilatation coefficient. Take an arbitrary vector m ^ 0. Let 

I: X = x{t), t G [tojTi] be an arbitrary curve issuing from the point xq, let 
x{to) = Xq] m be the tangent vector to I at xq, and let Ir* ■= x{[to,t*]) for 
each T* e (to, n]. 

The length dilatation coefficient at the point xq in the direction m is com- 
puted with the use of (3.4): 

Ki{xo,m,a,(3) := im 77— p = im °. ^ , , = . ' ■ 

The vectors x'{to) and m are coUinear, and hence x'^{to) = cm for some c. 
Eventually, we obtain 

Ki{xo,m,a,P) ^ - — p^^. (3.7) 

The volume dilatation coefficient. Take an arbitrary Vxo G V . Let v: x = 
x{S,), ^ G oj C R^, be an arbitrary three-dimensional subdomain of V containing 
Xq, Xq = x{^o); suppose that CO* Cw contains ^0 and := x{w*). 

The volume dilatation coefficient at the point xq is computed with the use 
of (3.6): 

Ks{xo,aJ):= lim ^= lim IfM^^M^^ = ^B^^ . (3.8) 

The volume in V^^^^^-^ of the parallelepiped spanned by k{xo) is equal to 

\Hxo)\^(xo) = l^'^l' and similarly, \k{xo)\a(^^„) = 

We denote aJ^ =: A; then a ~ A'^ f3, and hence the a- and /3-coordinates 
of an arbitrary vector y are related by y" = A~-^y^. Therefore, \k\a = |A:"| = 
= • Substituting this into (3.8), we obtain 

Ks{xo,aJ) = \A\ = \a{xo)\^^,^y (3.9) 
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The area dila,tation coeffieient. Take an arbitrary Vxo G V. Let s: x = x{u), 
u — (ui, U2) € cr C R^, be an arbitrary oriented surface in V with unit positive 
/3-nornial n at the point xq, xq = x{uo); we assume that a* C a contains uq 
and := x{a*). The area dilatation coefficient at the point xq is computed 
with the use of (3.5): 



K2{xo,n,a,$) := Um<^.^„o(|s<^. |^/|s<^. 1^)= 




(3.10) 



From this formula, we conclude that the dilatation coefficient is same as if 

1. (3{x) = /3(xo) and d{x) = a{xo)„ i.e., both metrics are Euclidean. 

2. s lies in the plane 11 with positive unit /3-normal n at the point xq. 
But in this case the integrands in (3.10) are constants and there is no need 

to pass to the limit; i.e., 

K2{xo,n,aJ) = Pf^ Vs C H. (3.11) 

Consider the parallelepiped v spanned by the vectors m, I (parallel to 11), and 

n forming a nondegenerate right frame at xq. Wc denote the parallelogram 
spanned by m and Z by s and the positive unit a-normal to 11 by n, 

= \s\fi ■ \n\^ = \s\^, \v\a = |s|a • (n,n)a. 

Therefore, in view of (3.11) and (3.9), 

K2{xo, n, a, 13) = -j- ■ (n, n)e, = |d|^ • (n, n}^. (3.12) 

Let us find n. To this end, we transform the equation of the plane 11 in the 
/3-metric to the equation in the d-metric. 

The equation 11 in the /3-metric has the form 

{n,x)fj = 0, (3.13) 

where x — ,to =: x, (n,a;)^ = {n^,x^). x" = A^^x^ , where A = oP . Therefore, 

{n,x)^ = (n^,y4x") = {A^n^,x'^). If we now introduce the vector TV with 
d-coordinates 

iV" := A^ni^, (3.14) 
then {n,x)p = {N",x") = {N,x)a. 
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By substituting this into (3.13) and by normalizing N in the d-norm, we 
obtain the equation of 11 in the d-metric: 

{N/\N%x)^=0. 

Hence n = iV/|A''"|, and using (3.14), we obtain 



(3.15) 



Now, using the expression of the d-coordinates via the /3-coordinates and 
(3.15), we have 

{n,n)a = (n",n") = {A-^n^ , A'^n^)/\A'^n^\ 

= H?/|aV| = i/|aV|. 

By substituting this into (3.12), we obtain 

K2{xo, n, a, /3) = , ' , (3.16) 

where A = a{xo)^'^''°^ ■ 



4 The strain frame. The relative stresses 

Material objects (curves, surfaces, three-dimensional regions) are deformed 
with respect to the initial state as a result of the displacement x = x{x). One 
can characterize the resulting strain by comparing the measures (lengths, areas, 
volumes) of these objects, as well as the angles between material curves, at time 
t and the initial time to. Therefore, we characterize all above-mentioned objects 
in the strained state, i.e., at time t, both by their geometric measure (length, 
area, etc.) at time t and by their initial measure at time to- 
Clear ly, the specification of a;(x) or x'(x) completely determines the terminal 
strain of the medium. In turn, the operator x'{x) =: A{x) at each point is 
determined by the specification of an arbitrary nondegenerate frame a(x) and 
the frame 

a{x) := (^(x)ai(x),^(x)a2(x),.A(x)a3(x))^ =: ^(x)a(x); (4.1) 

one even does not need to require that the function a(x) is continuous. 

However, if the frames a(x) are orthonormal, then, to specify the strain of 
the medium at time t, it suffices to define the frame d{x) alone at each point 
X G Vt- (Note that a{x) = a(x) for x = x{x); i.e., a{x) and a(x) specify the 
same function of a material point expressed via the Eulerian and Lagrangian 
coordinates, respectively.) 
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4-2. Let us prove this. The computation of the geometric measure of the 
deformed (i.e., given at time t) material objects does not encounter difficulties. 
Therefore, assuming that the frame field a{x) is given, we focus our attention 
on the computation of their initial measure. 

Let 1, Ij (i = 1,2), s, and v be oriented material curves, a surface, and a 
three-dimensional region given at the initial time in V; let I := a;(l), k := x{li), 
s :— x{s), and ?; := x{v) be the same objects at the terminal time; we suppose 
that li and I2 meet at the point xq; finally, xq = x(x.o); x = x(t); x = Xj(r), 
T G [ti,T2]; X = x(u), w(wi,U2) e fT C R^, and x = x(^), ^(^1,6,6) e w C 
are parametrizations of 1, 1^, s, and v. 

To denote the geometric measure of (length, area, or volume) of 1-3-dimensional 
objects, we shall use vertical bars. The initial measure of such objects given in 
the strained state will be denoted in the same way but with the subscript t. For 
example, |1| is the geometric length of 1, \v\ is the geometric volume of v, and \s\t 
is the geometric area of s (i.e., |s|). An geometric angle between vectors n and 
m will be denoted by {npih), and the initial angle between the same vectors 
issuing from a point x will be denoted by {npm)x,t- The adjective "geometric" 
for a measure will usually be omitted. 

First, let us compute the initial angle between the curves li at the point xq. 
For the parametrization of li one can take x = x{xi{T)) =: ^^(t). Let Tq be the 
parameter value corresponding to xq on the curves Ij; then it also corresponds 

to Xq on li. 

The derivatives a;-(ro) =: rrii are the positive tangent vectors to the curves 
k at the point xo, and x-(ro) =: rtij is the tangent vector to k at the point xq; 
here rrii and mj are related by 

rrii = Xi{To) = a;^(xo)x-(To) = ^(xo)mi. (4.2) 

We need the following lemma. 

Lemma. If A is a linear invertihle operator on E, & is a nondegenerate 
frame, and a = Aa., then for any at, (3 € E and an arbitrary frame b one has 

= (Aaf (4.3) 

< a,(3 >s,=< Aa,Af3 >a, |a|a, = |^a|a (4.4) 

|6|a = \Ab\a (4.5) 

Proof We have a = afs^. Therefore, Ace = afAai = afaf, i.e., {Aa.)"' = a**. 
Assertions (4.4) readily follow from (4.3) by virtue of definitions (2.2) and (2.3). 
Further, using the Gram determinant, we have 

|6|a = I (6^,b,)a|'/' = I {Abi,Abj)^\'^^ = \Ab\a, 

and the proof of the lemma is complete. 

Using (2.1)-(2.3), we compute the cosine of the initial angle between the 
curves k at the point xq- 

/ X / — ^ (mi, 012) (mi,m2)a(xo) 

cos(mi,m2)a:o,t = cos(mi,m2) = - 



|mi|-|m2| |mi|g,(xo) • |m2|a(xo) 
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Thus, using (4.4) and (4.2), we obtain 

cos(mi, m^)xa,t = 1 i i— V — • (4.6) 

l"*l|d(xo) ■ l"^2|a(xo) 

It is easy to express the initial length of the curve / via its equation 

= |1| = r |x'(r)| dr = r |x'(t)U(x(.)) dr. 
Since x'{r) = a;x(x(r))x'(r) = A{x{t))x' (r)^ we obtain, in view of (4.4), 

\l\t= r W{T)\^ixir))dT. (4.7) 
Jti 

In a similar way, we find the expressions for the initial area and volume: 



l^lf = |S| = / \/«i,<i> • «2'<2) - {<1'<2)^ du. 



Let us replace the scalar product (•, •) by (•, •)a(x)- (Recall that the frames a(x) 
are orthonornial.) Then we apply (4.4) and take into account the fact that 
Ax'^. to obtain 



The coordinate frame of the curvilinear coordinates x = x(^) at the point k(x) := 
f dx dx dx\^ dx r if / m ,^ r,- 

[Wi'dC^'d^J d^' ' ^ ^ I""'* " '"'I " Ja;|k(x)|d6 Smce the 

frames a(x) are orthonormal, it follows that E = ^^^(x)! and, using (4.5), we 
obtain |k(x)| = |k(x)|a,(^) = |^(x)k(x)|a(^). Since 



^(x)k(x) = 4(x)^ = (4(x)x^^,4(x)x^^,4(x)x^3)^ =: k{x) 

is a coordinate frame of the curvilinear coordinates x = x{^) at the point x, we 
have 

\v\t= I |fc(a;)U(,)de (4.9) 

Formulas (4.6)-(4.9) show that, given the frame field a{x), one can compute 

the initial measure of 1-3-dimensional material objects and the initial angle 
between the material curves from the equations of these objects in the strained 
state (i.e., at time t). It follows that by specifying the field a{x) one completely 
determines the strain of the medium at time t. 

On the other hand, Eqs. (4.6)-(4.9) show (see Section 3.1) that the origi- 
nal metric for deformed material objects is the metric generated by the frame 
field a{x), i.e., the metric obtained in the Riemannian afiine space Vt^a by the 
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replacement of the scalar product (•, •) common for all points in Vt by the field 
of scalar products (•, •)a(a;)) individual for each point x € Vt- 

4-3. We take a point x &Vt and a vector m ^ 0. Let I be a material curve 
at time t for which m is the tangent vector at the point x. We consider the 
length dilatation coefficient for this curve at time t, equal to the ratio |^|/|^|f of 
its length at time t to its length in the initial state. 

It turn out (see Section 3.2) that the limit of this coefficient as the curve 
shrinks along itself to the point x, is independent of the choice of I and depends 
only on x and m. We refer to this limit as the length dilatation coefficient at 
the point x at time t (i.e., at the point {x,t)) in the direction m and denote it 
by Ki{x,m). 

The area dilatation coefficient K2{x, n) at a point (a;, t) of a material surface 
with unit normal n at the point x at time t and the volume dilatation coefficient 
Ks^x) at the point {x,t) are introduced in a similar way. 

Since for each orthonormal frame e the scalar product •) coincides with 
(•, •)e, it follows that the geometric measure on Vt coincides with the measure 
on Vt^e for any orthonormal frame field e(x) on V*. Furthermore, the initial 
metric coincides with the metric in Vt^a- Therefore, Ki is the measure dilatation 
coefficient for the passage from the a- metric to the e- metric (i.e., the geometric 
metric) in I4. Hence, according to formulas (3.7), (3.16), and (3.9) with d = a 
and $ = e, we have 

i^i(x,m) = ^, K2{x,n) = j^^, K^{x) = \d\, (4.10) 

where A = d^. 

Formula (4.6), together with (4.10), shows that a{x) completely determines 
the strain of the medium at the point {x,t). 

Therefore, d{x) will be called the strain frame. Let us give a complete defi- 
nition taking into accoimt (4.1) and using derivatives with respect to a vector. 

The strain frame of the medium at a point x at time t is defined as 

- ^ - = (^'(^)- ^'w-)^ 

= a;'(x)a, 

(4.11) 

where a is an arbitrary orthonormal frame. 

The physical meaning of the strain frame d is clear from the very definition 
of the derivative with respect to a vector in conjunction with (4.10) (which is 
in generally not necessary) for m = tti. 

The vectors of the frame d are the tangents at the point x to three strained 
material fibrils issuing from x and mutually orthogonal at the initial time. More- 
over, \ai \ is the dilatation coefficient of the material at x at time t in the direction 

Restating the last paragraph in Section 4.2, we can say that the specification 
of the strain frame field permits one to introduce a Riemannian metric on Vt 
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such that the measurement of angles between deformed material curves and of 
lengths, areas, and volumes of deformed 1-3-dimensional material objects in 
this metric gives the values of the respective characteristics for the objects in 
question prior to strain without resorting to their initial state. 

The strain frame is not unique. By taking another orthonormal frame a' in 
the undeformed medium, we obtain another strain frame a' := Aa.' by (4.11). 
Obviously, the frames d and d' characterize the same strain. 

Note that the values of the expressions (4.6) and (4.10) do not change if the 
frame a is replaced by an equivalent frame. 

4-4- Wc take a point x in the deformed medium and a unit vector n. The 
stress per unit geometric area at the point x at time t on the oriented plane with 
positive normal n will be called the geometric stress (the Cauchy stress [12], or 
the true stress [13]) and will be denoted by (Tn{x). The stress per unit initial 
area will be referred to as relative (to the initial area). 

The geometric measure on Vt coincides with the measure on Vt,e for any 
orthonormal frame field e{x) on Vt, and the initial measure, as was already 
mentioned, coincides with the measure on Vt^a- 

Generalizing, for an arbitrary nondegenerate frame field b{x) defined on Vt, 
we introduce the stress Un{x,b) per unit area in the Riemannian space V^f^, 

which will be called the relative stress (relative to b{x)). The symbol b is omitted 
if it is clear what frame is meant. 

The geometric stress and the stress relative to the initial area are the special 
cases of Un{x,b) for b{x) = e{x) and b{x) = d{x), respectively. 

4-5. We introduce the following notation. Let bi be one of the vectors in the 
frame b. By ni{b) we denote the unit normal to the plane spanned by the other 
two vectors of b; we assume that ni{b) points to the same side of this plane as 
bi. Next, Ui{x,b) := u^,^jj^^{x,b) are the relative stresses on the faces of the 

frame b{x) at the point {x, t); n is the unit positive normal to a plane 11 passing 
throTigh x: n is the unit positive normal to the same plane in the metric V^ j (it 

will be called the 6-normal); (ni,n2,n3)^ := n''. 
One can prove the generalized Cauchy relation 

Un{x, b) = Ui{x, b)ni (4-12) 

for the relative stresses in completely the same way as the classical Cauchy 
relation between the Cauchy (geometric) stresses on the plane 11 and the faces 
of e for an orthonormal frame e. 
Introducing the frame 

u{x,b) := {u-i{x,b), U2{x,b), 113(3;, 6)) (4.13) 

of relative stresses on the faces of b{x), we can rewrite the dependence (4.12) in 
the form 

Un{x,b) ^u{x,bfn''. (4.14) 

4-6. Let us establish the connection between the relative stress frames for 
different initial frame fields. We take a point x G Vt and two arbitrary nonde- 
generate frames a and b at x. (The values of the frames at any other points are 
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irrelevant.) Suppose that u{x,a) is known. Let us compute the stress Ui(x,b) 
on the ith face of b. We denote the unit positive normal, the d-normal, and 
the 6-normal to this face by m, n, and n, respectively; next, A := a^, B := b", 
and i is the column with unit in the ith position and zeros in the remaining 
positions. Then n = bi and hence 

= i; (4.15) 

since d = A'^b, it follows that b = A~'^d and hence 

B = A-^. (4.16) 

By formulas (3.15) and (3.16), the area dilatation coefficient for the transi- 
tion from the a-measure to the 6-measure is equal to 

K2{x,n) = -^ and n» = (4.17) 



Using (4.17), (4.14), and (4.15), we obtain 

-u{x, d) n' 



Ui(x, b) = K^^ix, n)um(x, d)= — — — u(x, d)"^n"= 
= \d\^'^u{x,d)'^A'^i = \d\^^u{x, d)'^ {an, ai2, ais)'^ = \d\^^Uj{x,d)ai 



and since, according to (4.16), the 6-volume |d|j, of the parallelepiped spanned 
by d is equal to the determinant \A\ = \B\~^, we have 

u{x,b) = \B\B-^u{x,d). (4.18) 

This is the desired relation between the relative stresses. 

Let us clarify the relationship between the matrice Ua of decomposition 
of u{x,d) with respect to the frame d and matrice 11^ of decomposition of 
u,(x, b) with respect to the frame b. Successively substituting u,(x, d) = U]^d 
and d = B^'^b into (4.18), we have u{x,b) = \B\B'^U^ B^'^b and hence 
{/j = \B\B~^UaB~^ . Thus the matrix of decomposition of the frame of rel- 
ative stresses on the faces of the frame d with respect to d changes under the 
replacement of d by 6 as the matrix of some contravariant pseudotensor of rank 2 
and weight 1. It will be called the relative stress pseudotensor at the point x. 



5 The equilibrium equations with respect to 
an arbitrary frame field 

Suppose that a nondegenerate frame field b{x) of class is given in V and an 
orthonormal frame e is fixed. Let F{x) and ^(x, &) be the mass force density 
at the point x per unit volume (geometric density) and per unit relative b- 
volume (relative density), respectively, and let (Ti{x) •.= Ui{x,e) be the Cauchy 
(geometric) stress at x on the ith face of e. 
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Since the volume dilatation coefficient (3.9) for the transition from the b- 
volume to the geometric volume is equal to we have *(a;, b) = F{x)\b\. 

The force and moment equilibrium equations for geometric stresses (e.g., see 
[13]) read 

^^+F{x) = 0, S = S^, (5.1) 

where S := and the Xi are the coordinates of x in the Cartesian coordinate 
system with coordinate frame e. By (4.18), we have 

a = \b\-^Bu, (5.2) 

where B = b^ and hence (Ti{x) = \b{x)\~^bij{x)uj{x,b). Therefore, 

The first term in (5.3) is equal to 

\b\~^Uj divbj. (5.4) 

Let Sij be the Kronecker delta, and let b^^ be the entries of the matrix B~^. 
Using the linearity of the derivative with respect to a vector in the vector and 
the fact that b = B'^e, i.e., 

e = B-^b (5.5) 
and hence = b'^^{x)bk{x), we see that the second term in (5.3) is equal to 

By substituting (5.4) into (5.3) and the resulting relation into the first equation 
in (5.1), we obtain the force equilibrium equation for the relative stresses: 

duj { d\n\h\\ , , 

Let us compute the matrix =: U. By substituting the expression (5.5) 
into the equation a = S-^e and then the resulting relation into (5.2), we have 
u = \b\B-'^Y,'^ B-'^b. We see that the matrix U is equal to \b\B-'^i:B-'^ and 
hence (sec (5.1)) is symmetric. 

Therefore, the moment equilibrium equation for the relative stresses is also 
reduced to the symmetry condition for a matrix, in this case, for the matrix U : 

U = U^. (5.7) 

Remark. The arbitrary nondegenerate frame field b{x) occurring in Eq. (5.6) 
need not be the coordinate frame field of any curvilinear coordinate system. 
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6 The elastic strain frame. The elastic state 
equation 

6.1. We momentarily assume that the medium in question is clastic, was in the 
natural (unstressed) state at the initial time, and was subjected only to elastic 
strains until time t. In this case, the initial metric of the deformed medium will 
be referred to as the natural metric; accordingly, the stress Un{x,a) per unit 
natural area will be called the natural stress (it is also known as the conditional 
stress vector [13] and the first Piola-Kirchhoff stress vector [12]) and will be 
denoted by tn{x). The density \l/(a;, a) of mass forces per unit natural volume 
will also be called the natural density and will be denoted by $(x). 

Now the strain frame defined by formula (4.11) completely characterizes the 
elastic strain of the medium at the point {x,t). It will be called the elastic strain 
frame. 

Since only natural stresses and densities are used in what follows, we usually 
omit the word "natural." 

The Cauchy relation (4.12) passes for the natural stresses into the well-known 
relation [12, 13] 

tn{xo) =ti{xo,a)ni, (6.1) 
where ti{xo, a) := t.^f^Q^^{xo), n is the unit a-normal at the point x to the plane 
orthogonal to n, and (ni,n2,n3)^ := n". 
Recall that, according to (4.17), 



A = o^ (6.2) 



for b = e. 

We assume that the medium in the natural state is homogeneous and isotropic 
and that the stress at an arbitrary point x of the deformed medium is completely 
determined by strain at x, i.e., by the elastic strain frame d{x). 

Using the frame i{x, a) of natural stresses at the point x on the faces of 
the frame d{x), we introduce the stress matrix t{x,d) := {tij{x,d)) := {""{Xjd). 
Under our assumptions about the medium, the strain frame d{x) determines the 
stress matrix by the following chain of relations [12]: d e'^{d) := {efj{d)) := 

{^{{ai,aj) — Sij)) (the Cauchy-Lagrange elastic strain matrix); /(e*^) := 

{fij{e'^)), where are functions of class characterizing the elastic properties 
of the medium; t{x,d{x)) = f{e^{d{x))). 

Thus the equations of the elastic state of the medium can be represented in 
the form 

t{x,d{x)) = f{i%d{x))), (6.3) 

or 

U{x,d{x)) = fj^{i''{d{x)))aj{x), i = 1,2,3. (6.4) 

Thus the matrix i{x,d), which is a special case of U (see Section 5), is 
symmetric, and the principal stress frame is coaxial with the principal strain 
frame. 
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Since in our case the functions t, i, and ti depend only on the frame a{x), 
the first argument x of these functions can be omitted in what foUows. 

6.2. Now consider an elastoplastic medium in the process of elastoplastic 
deformation. First, we assume that, at any time t of the process, the medium, 
or at least a sufRciently small neighborhood of any point of the medium, can 
be unloaded to a natural state by an elastic deformation. We treat this natural 
state and the state at time t as the initial undeformed state and the state 
obtained from it by an elastic deformation, respectively. Then (see 6.1) there 
exists an elastic strain frame at time t at each point of the medium. 

There also exists a similar frame field under more general assumptions. In 
general, even a small neighborhood of a point cannot be unloaded to an un- 
stressed state. There always remain residual stresses. However, one can assume 
that these stresses tend to zero as the neighborhood shrinks into the point, and 
so the unloaded states tend to an unstressed state. 

We take an arbitrary point x of the deformed medium at time t. Let v be 
a neighborhood of x. Let us unload the part of the medium contained in v and 
consider the strain frame dy{x,t) at x, corresponding to the displacement of 
the unloaded part to its position at time t. Suppose that there exists a limit 
limay{x,t) =: d{x,t) as w — > a; and that the field d{x,t) is continuous. 

Now consider an arbitrary material curve in the deformed medium. We cover 
it by finitely many balls of radius r and cut into finitely many partial curves 
each of which entirely lies in one of the balls. By unloading the part of the 
medium lying in a ball, we obtain the " unloaded" length of the corresponding 
partial curve. We define the "unloaded" length of the entire curve as the sum 
of unloaded lengths of its parts. 

One can show that the unloaded lengths of the curve tend as r ^ to a limit, 
which will be called the "natural" length of the material curve in question. The 
natural measures of material surfaces, volumes, and angles can be defined in a 
similar way. It turns out that the natural measure defined thus constructed for 
material objects coincides with their measure in the Riemannian space generated 
on Vt by the frame field d{x,t). 

Thus in both cases at any time t at each point of the medium there exists a 
frame determining how the measure of 1-3-dimensional material objects at time 
t has changed compared with their natural measure. It is this property of the 
frame that is important in what follows. 

Therefore, let us introduce a notion of elastic strain frame for elastoplastic 
media generalizing the corresponding notion for elastic media and based neither 
on any concept of actual motion of the medium nor on any assumptions related 
to full or partial unloading of parts of the medium. 

Suppose that at each time the medium has some " elastic structure" , which 
permits one to judge about the natural measure near each point of the medium. 
More precisely, we shall adopt the following assumption. 

Assumption. For each given time t, to each point x of the medium one can 
assign a nondegenerate frame d{x, t) (which will be called the elastic strain frame 
of the medium at the point {x^t)) such that the Riemannian metric generated 
in Vt by the field d{x,t) (the metric of the space Vt,a) is "natural." 
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Remark. 

1. The "naturality" of a measure, which is clear if full unloading is possible, 
is defined as certain relationships, postulated below, with the stress and the 
plastic strain. 

2. In a polycrystal material, crystal lattices can serve as a prototype of an 
"elastic structure." 

The clastic strain frame field on Vt is not unique. Indeed, if d{x,t) is an 
elastic strain frame field, then for a frame field d'{x,t) to be an elastic strain 
frame field it is necessary and sufficiently that these fields be equivalent, i.e., 
that the frames a{x,t) and d'{x, t) be equivalent for all .x £ Vj. 

Indeed, if the fields d and d' are equivalent, then these frames generate equal 
scalar products at each point. Consequently, the Riemannian metrics generated 
by the fields d and d' coincide: i.e., d! also generates a natural metric and hence 
is an elastic strain frame field. The necessity can also be easily justified. 

The following lemma will be often used in the sequel. 

Lemma. For nondegenerate frames d and d' to be equivalent, it is necessary 
and sufficiently that there be an orthogonal 3x3 matrix Q such that 

d' = Qd. (6.5) 

Indeed, let the frames be equivalent: Va, /3 {cx.,l3)a' — {a.,l3)a. 

We denote Q := (d'°)^; then d' = Qd. Since each of the frames is or- 
thogonal with respect to the scalar product generated by it, it follows that 
Sij = {a[,a'j)a' {a[,a'j)a = {qikak,qjiai)a = qikq]i{ak,a,i)a = qikqjiSki for 
any i and j, i.e., qikQjk = ^ij and hence Q is orthogonal. The sufficiency can be 
proved even easier. 



7 The plastic strain tensor 

To characterize plastic strain, we need additional information on the change in 

the measure of geometric objects midcrgoing deformation. 

7.1. The measure dilatation coefficients for the mapping of Ya into Vg. Here 
the situation is the same as in Sections 4.2-4.3; Ya plays the role of the domain 
of the initial position of the medium, plays the role of the domain Vt of the 
final position of the medium, d{x) is used instead of a(x), and the /3-measure is 
used instead of the geometric measure in Vt . We introduce the frame 

w N 9x{x) 

which plays the role of d{x) in 4.2-4.3. Arguing as in 4.2-4.3, we arrive at the 
conclusion that the d - measure of objects in V coincides with the d-measure of 
their preimages in V. But then the following assertion holds. 

Assertion. The measure dilatation coefficients comparing the d-measure of 
1-3- dimensional objects in V with the ^-measure of their images in V are given 
by formulas (3.7), (3.9), and (3.16) with a replaced by d. 
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For the a-mcasure of the initial angles, using the reasoning and notation of 
Section 4.2, we have 

cos(mi,m2)x,a = 1 f . — V — , (7.1) 

where (mi7in2)x,a is the angle between mi and m2 at the point x in the a- 
metric. 

7.2. In what follows, wc consider a medium that moves during the time 
interval [tb,te\ =■ T and is elastoplastic (in the sense of Section 6.2) at each 
time t G T. All earlier-introduced assumptions and the main notation are 
preserved; wc only introduce an additional variable, time {t), which occurs as 
an argument of the functions specifying the law of motion of the medium, the 
elastic strain frame, the velocities v{x, t) of points of the medium, etc. 

Let W := {{x,t)\t <= T, x E Vt}. Let a{x,t) be an clastic strain frame field 
of the medium of the class C^{W). We take an arbitrary to GT and denote the 
law of motion of the medium starting from time to by 

x = Xto{x,t). (7.2) 

To each point x in space, this law assigns the spatial position x at time t G T 
of the material point residing at x at time to- We shall assume that xtg G 

C'{Vt, X T). 

The comparison of natural measures of material objects at time t and to 
characterizes the plastic strain of the medium on the interval [to, t]. 

Just as the frame field dx/dh on Vj, where a defines the geometric measure 
at the initial time (since E = Eg^), completely characterizes the change in the 
geometric measure of material objects (see Section 4.3), so the frame field 

^(^'*)^=l|!;^' (- = --(x,0), (7.3) 

where a{x,to) defines the natural metric on Vtg at time to, completely charac- 
terizes the change in the natural measure of material objects in time [to, t], i.e., 
the plastic strain on this time interval (the elastic strain frame field d{x, t) on 
Vt and VtQ being known). 

Indeed, we take Vxo G Vi„ and xo := xi^ (xo, t) and momentarily set a(xo, to) =: 
do =: (01,02,03)^, a{xo,t) =: d, and b{xo,t) =: b. According to the Assertion 
in Section 7.1 (with a as and 6 as d), the natural length dilatation coefficient 
at a point {xo, t) in an arbitrary direction m is equal to 

K,{xo,t)=^-^. (7.4) 

Moreover, for two arbitrary nonzero vectors mi and m2 at xq, one can take 
oriented material curves for which these vectors are the positive tangent vectors 
at the point xo at time t and introduce the positive tangent vectors rrij to these 
curves at time to at the point xq. 
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Now, using the properties of the elastic strain frame d (sec the assumption 
in Section 6, (4.6), and (7.1) with d = ao), for the material angle [11x^1712) xo,t 
and the corresponding (iri]7m2)xo,to we have 

. (rui, 7712)0 

cos(mi,Tn2)^^_t = 



(7.5) 

mi,m2)d„ _ {mi,m2)i 



\m1\a\m2\a 

cos (mi'm2)^ t = , , , , — , 

^ ''^o,to |miUJm2U„ |mi|^|m2|j 

Thus the specification of b for known d permits one to compute the natural 
length dilatation coefficient at xq at time t in any direction m (see (7.4)) and, 
for any material curves issuing from xo at time t, compute the cosines of the 
natural angles between them at the point xq at time t and at the point xq at 
time to (see (7.5)). 

In particular, for the vectors 6j of the frame b, we have b\ = i and \bi\i = 1. 
Hence (see (7.4) and (7.5)) \bi\a is the material length dilatation coefficient along 

bi at the point xo at time t; moreover, cos{bi,bj)xo,t = {bi,bj)al{\bi\a\bj\a), 
while the cosine of the material angle between the corresponding vectors a° and 
aPj of the frame ao at xo at time to (see (7.5)) is zero for i ^ j. 

We see that the frame b plays the same role in determining the plastic strain 
of the medium in time [to,t] as d plays in determining the elastic strain. 

By analogy with the Cauchy-Lagrange elastic strain matrix, we introduce 
the plastic strain matrix at a point xq from time to- 

iP:=(li{b,,bjU-Sij))=:^. (7.6) 

By introducing the matrix B := 6*^ and by taking into account the definition of 
the scalar product (•, •)„, we obtain 
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The physical meaning of is the same as that of with the only difference 
that geometric lengths and angles arc replaced by natural ones. Let us study 
how sP changes if we replace the elastic strain frame ao at the point (xq, to) by 
the equivalent frame Oq. It suffices keep track of the change in 7. We denote 
Qo := {{a'of'V; then 

do = Qodo- (7.7) 

For the objects b' , B' , (e^)', and 7' corresponding to do to belong to the 
class C^, we replace the field a{x,t) by an equivalent field d'{x,t) such that 

d'{x,t) = Q{x,t)d(x,t), (7.8) 

where Q{x,t) is a field of orthogonal matrixes (see the lemma in Section 6) 
and Q{xo,to) = Qo- Now, in view of (7.3), (7.7) and (7.8), we obtain 



9d'(xo,to) 



-d'o = Qob = QoB' d = QoB' d'. 
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It follows that 

B' = QBQ'^. (7.9) 

Therefore, 7' = B'^B' = QqB^ BQq = QotQo i i-6-> 7 is the matrix of some 

covariant rank two tensor 7 in the frame do. This tensor is defined at (xo,fo) 

for elastic strain frames and depends on time. 

Hence the same it true for and for the tensor of plastic strain at the 

point xo from time to- 
la general, the present paper deals only with rank two tensors, and since 

all of them are associated only with elastic strain frames, it follows that the 

matrices of transformation from one frame to another are orthogonal, so that 

the tensors are covariant. 

Remark. The tensor is different from the tensor, traditionally denoted 

by the same letter, defined as the difference between the full and elastic strain 

tensors or reckoned from the initial time. 

Let us introduce the plastic strain rate tensor ^(xo,to) at the point (xo,io)- 

By definition, it is equal to the total i-derivative at time to of the plastic strain 

tensor at the point xq reckoned from time to, 

C(xo,to) = (d£7di), . (7.10) 



8 The strain consistency equation 

On the mass point trajectory passing through xq at time to, we have 

^^*g^°'^^ "(^0'^o) = iB{x,t))'^a{x,t), x = Xio(xo,t) 

according to (7.3). By taking the total ^-derivative of this identity at i = to and 
by changing the order of differentiation on the left-hand side, we obtain 



dv(xo,to) 
9a(xo,to) 



d(xo,to) + -B^(xo,to) 



-a{x,t) 



Taking into account the fact that B'^{xo,to) = I and passing to contracted 
notation, we obtain the equation 



dv 

da 



fdB^ 



\ dt 



dt) 



(8.1) 



satisfied at the point (xo,to) by any clastic strain frame field d and by the 
corresponding B (since we deal with a specific motion of the medium, the field 
V is fixed). 

By taking the total t-derivative of (7.9), for t = tg we obtain the law describ- 
ing how the matrix dB/dt\^^ changes under the passage to an arbitrary elastic 
strain frame field d'{x,t) (see (7.8)): 



dB' 



(8.2) 
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For this law to be a tensor law on W, it is necessary and sufficient that 



^9(^ = oyM)ew, (8.3) 

at 

i.e. that Q{x,t) be constant along the space-time trajectories of the points of 
the medium. 

We say that elastic strain frame fields a and a' related by formula (7.8) with 
Q satisfying (8.3) arc kinematically equivalent. 

Thus the set of elastic strain frame fields splits into classes of kinematically 
equivalent fields, and the matrix dB/dt\^^ varies according to the tensor law on 

each of the classes. 

These classes arc sufficiently ample; namely, the following assertion is valid: 
if a{x, t) is an elastic strain frame field of the medium in W , then for each tg and 
each clastic strain frame field aQ{x) given at time tg on Vta there exists a field 
a*{x, t) kinematically equivalent to a{x, t) and coinciding with aa{x) at time ig- 

Indeed, let Qo{x) be the field of orthogonal matrices such that ag(x) = 
(5g(x)a(x, tg). We introduce a field of orthogonal matrices on W by setting 
Q{x,to) = Qo{x) and by requiring that Q{x,t) docs not vary along the tra- 
jectories of the points of the medium. Let d*(x,t) := Q(x,t)a{x,t)\ then first, 
a*{x,tQ) := Q{x,to)a{x,to) = ao{x) and second, the field a* is kinematically 
equivalent to a{x,t), since dQ{x,t)/dt = 0. 

Assumption. For the medium, there exists an elastic strain frame field of the 
class C^{W) for which the matrix dB/dt\^^ is symmetric on W. We say that 
this field is kinematically constrained. 

Remark 8.1. It is shown in the paper mentioned second in the footnote 
in Section 1 that if the law of motion of the medium is of the class , then 
the existence of a kinematically constrained field automatically follows from the 
existence of at least one clastic strain frame field of the class C^. 

It turns out that the set of kinematically constrained elastic strain frame 
fields is a class of kinematically equivalent elastic strain frame fields. 

Indeed, if a is a kinematically constrained elastic strain frame field and 
a' is a kinematically equivalent field, then {dB'/dt)\^^ = Qo{dB /dt)^^^QQ by 
(8.2), and hence this matrix is symmetric (since so is dB/dt\^^); thus a' is 
also kinematically constrained. On the other hand, if elastic strain frame 
fields a and a! arc kinematically constrained, then the matrices dB/dt\^^ and 

dB' /dt\^^ are symmetric, and hence so is (^dQ / dt\^^ Qq (see (8.2)). Conse- 
quently, {dQ/dt)\^^QQ = Qo{dQ^ /dt)\^^, and hence 

di|,„^° 2 Uilto ° UoJ 2 dt Uo 

since QQ^ = I. Thus dQ/dt\^^ = 0; i.e., a and a! are kinematically equivalent. 

The class of kinematically constrained elastic strain frame fields plays a 
special role. Let d be a kinematically constrained elastic strain frame field, and 
let (xo,fo) be an arbitrary point in W. Taking into account the symmetry of 
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dB/dt\^^^, the relation B\^^ = I, and formulas (7.10) and (7.6), we obtain the 
following expression for the matrix C(xo, ^o) in the frame ao: 



Cao(xo,*o) = 



dt 



dt 



fdB^ 



V dt 



B 



I to 



B' 



,dB 
'dt 



ItQ 

dB^ 



I to 



Substituting this into (8.1), we obtain the strain rate consistency equation 



dv 

da 



C.aa + 



da 

dt' 



(8.4) 



which is satisfied in W by each kinematically constrained elastic strain frame 
field. 

The physical meaning of Eq. (8.4) becomes clear if we remember how (8.1) 
was obtained: the material rate dv{x,t)/da{x,t) due to the motion of the 
medium of the frame coinciding with d{x,t) at the point {x,t) at time t is 
equal to the material rate da{x,t) / dt of the elastic strain frame plus a correc- 
tion term depending on the matrix of the plastic strain rate tensor in the frame 
a. 

Rem,ark 8.2. For "small" clastic strains, more precisely, under the assump- 
tion that there exists an orthonormal frame e and a number 5 « \ such that 
\ai — Bil and \dai/dxj\, \dai/dt\, \dv/dxi\ < S, it follows from (8.1), neglecting 
quantities of the order of 5^, that 



dt 



Sto {xo,to) 



^ e'tg{xo,to) + ^e''{xo,to), V(a;o,io) e W, 



dt 



where the full and plastic strains are reckoned from time to. It is this equa- 
tion that essentially substitutes the traditional equation de = de^" + de^. 
(Although this substitution is not justified, it leads to correct results in our 
situation.) 



9 The elastic and plastic constitutive equations 
for 

elastoplastic media " homogeneous and isotropic 
in the natural state" 

We understand the media specified in the title as media possessing the following 
properties. 

Assumption 9.1. Under elastoplastic deformation, the elastic strain frames 
and the natural stresses are related by the same elastic constitutive equation 
(6.3)=(6.4) as in the case of elastic media homogeneous and isotropic in the 
natural state. 
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Let a* be the principal clastic strain frame at a point {x,t) G W: then, 
by Assumption 9.1, the principal frame i{d*) = {ti{a*),t2{a.*),t3{a*))'^ of nat- 
ural stresses on the faces of a* and the frame a* are coaxial; here ti{a*) are 
the principal natural stresses. The only possible nonzero coordinate ti{a*) 
in the frame d* is tu{a*) =: ti{a*) will also be referred to as the principal 
stress, which, in view of our notation, will not lead to a misunderstanding; 
{h{a*),t2{d*),h{a*)) =:t(d*). 

Assumption 9.2. If d* is the principal elastic strain frame of the medium at a 
point {x,t), then the tensor C(^)*) is coaxial to d* for all {x,t) G W; moreover, 
there exist known fimctions pi : — > R (i = 1,2,3), describing the plastic 
properties of the medium such that 



Ca* (a;, t) 



Pi(t(d*)) 
P2(t(d*)) 
P3(t(d*)) 



=: P(t(d*)). 



(9.1) 



Remark. 1. Thus the medium in question is so far "perfectly plastic" in 
the sense that the plastic strain rates are independent of the strain history and 
depend only on the principal stresses. 

2. The assumption that Ca* depends only on the principal stresses and is 
independent of the stress increments seems to us to be quite natural for "slow" 
deformations. 

3. Using the technique applied in the study of the yield locus of an isotropic 
material (e.g., see [14]), one can readily show that 

Pi{t2MM) =P2{h,ti,t2) =Pi{ti,t2,ts) =:p(t) Vt. 

For an arbitrary elastic strain frame d at a point {x,t), by i? we denote the 

orthogonal matrix such that d = Rci* . Then Ca = RCa'R'^, and the plastic 
constitutive equation of the medium (9.1) can be represented in the form 



Caix,t) = RP{t{d*))R'^ y{x,t) G W. 



(9.2) 



10 The system of elastoplastic strain equations 
in the Eulerian variables 

Let us summarize our study. We consider an elastoplastic medium satisfying the 
assumptions in Sections 6 and 8 and moving in a space-time domain W. The 
natural density ^{x,t) of applied mass forces of the class C{W) and the den- 
sity po of the undeformed medium are known. The unknowns are the velocity 
field v{x,t) of the points of the medium and any of the kinematically equiva- 
lent kinematically constrained elastic strain frame fields a{x, t) of the medium, 
where d and v are related by two equations, namely, the strain consistency 
equation (8.4) (or (8.1) is d not kinematically constrained) and the equation of 
motion obtained from the force equilibrium equation (5.6) by the substitution 
of a{x,t) for b{x). Here the relative stresses Ui{x) become the natural stresses 
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ti{x,t), and the mass forces include also inertial forces; i.e., ^(x) is replaced 
by *(a;,i) - poidv/dt). 

These equations are valid for any media with an "elastic structure." Along 
with the unknowns d and v, they also contain the natural stress frame i and 
a matrix describing the plastic strain; therefore, to obtain a closed system, one 
should supplement these equations with equations describing the elastic and 
plastic properties of the material. 

For media homogeneous and isotropic in the undeformed state, i.e., satisfying 
Assumptions 9.1 and 9.2, the supplementing equations are the elastic and plastic 
constitutive equations (6.4) and (9.2). 

Thus the system of dynamic equations of elastoplastic strain in this case has 
the form 



dv dti / 91n|a| 

Po-TV = ^ ^ 

dt oai \ oai 

^ = ^-RPR^d, \f{x,t)eW 
dt oa 



(10.1) 



U = fj,{e%a))aj. (10.2) 

The total derivatives in (10.1) are expressed via v and d as follows: 

dv dv dv da dd dd 
'dt^lhi^'dt' dt^d^^'dt' 



(10.3) 



System (10.1) (in view of (10.2)) is a system of a vector equation and a frame 
equation for the vector function v and the frame function d. In vector form, it 
consists of four equations for v and {i = 1,2,3). 

The equations are represented in coordinate-free form. After the introduc- 
tion of some coordinate system and the replacement of derivatives along vectors 
by their expressions via partial derivative with respect to coordinates, one ob- 
tains a system of quasilinear first-order partial differential equations. 

Equations (10.1) are solved for the time derivatives of the desired functions, 
which makes the system suitable for numerical time-marching solution. 

If the functions Pi describing the plastic properties of the medium depend 
not only on the principal stresses but also on parameters determined by the 
strain history, then one should supplement the system with equations relating 
these parameters to the desired functions. 

In the absence of plastic strains {pi — 0), system (10.1) passes into system 
(41) of elasticity equations!! 
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